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Abstract 

We obtain general analytical forms for the solutions of the one-loop renormal- 
ization group equations in the top/bottom/r sector of the MSSM. These solutions 
are valid for any value of tan (5 as well as any non-universal initial conditions for the 
soft SUSY breaking parameters and non-unification of the Yukawa couplings. We 
establish analytically a generic screening effect of non-universality, in the vicinity 
of the infrared quasi fixed point, which allows to determine sector-wise a hierarchy 
of sensitivity to initial conditions. We give also various numerical illustrations of 
this effect away from the quasi fixed point and assess the sensitivity of the Higgs 
and sfermion spectra to the non-universality of the various soft breaking sectors. 
As a by-product, a typical anomaly-mediated non-universality of the gaugino sector 
would have marginal influence on the scalar spectrum. 
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1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) |l|] has been intensively studied 
recently with the emphasis on prediction of particle spectrum. It crucially depends 
on the mechanism of SUSY breaking and the one which is commonly accepted introduces 
the so-called soft terms at a high energy scale || ||. These soft terms are running then 
down to low energy according to the well known RG equations starting from some initial 
values. In the minimal version the soft terms obey the universality hypothesis which 
leaves one with a set of 5 independent parameters (before radiative electroweak symmetry 
breaking) : m , m^, A , /x and tan/3. However, recently there appeared some interest 
in new contributions to SUSY breaking patterns which result in non-universal boundary 
conditions for the soft terms ||. To investigate the influence of non-universality on the 
particle spectrum it is very useful to have analytical solutions to the RG equations for 
the evolution of the soft terms. Numerical analysis though straightforward is rather 
complicated due to a large number of free parameters while an analytical solution allows 
one to see which of these parameters is important and which is not. 

While with a single Yukawa coupling the analytical solution to the one-loop RG equa- 
tions has been known for long, for increasing number of Yukawa couplings it has been 
obtained quite recently || in the form that allows iterative presentation. On the other 
hand, it has been shown |7J that when one knows the solution to RG equations for the 
couplings of a rigid theory, one can obtain those for the soft terms by a usual Taylor 
expansion over the Grassmannian parameters. Thus one can apply the Grassmannian 
expansion to the general solutions of || to get the analytical solution for the soft terms 
with arbitrary initial conditions in an iterative form and explore their dependence on the 
latter. 

In the present paper we perform this analysis in the MSSM with three Yukawa cou- 
plings (Y t , Yb and Y T ) and arbitrary initial conditions for the soft terms in the one-loop 
approximation. We examine the dependence of the solutions on the initial conditions 
and show that if the Yukawa couplings are large enough (typically bigger than a s at the 
GUT scale) the solutions exhibit a quasi-fixed point behaviour. This means that in such 
a regime some of the initial conditions are completely washed out, and become actually 
inessential at low energy. The main role is played by the gaugino masses out of which the 
gluino mass is by far dominating. 



2 Analytic Solution to the RG Equations for Yukawa 
Couplings 

Though the RG equations for the Yukawa couplings do not have explicit analytic solution, 
they can be solved iteratively as it has been demonstrated in Ref. 0. Using the notation 

2 2 

on = -^t, i = 1, 2, 3; Y k = k = t,b,r 

16tt 2 ' ' ' ' 16tt 2 ' ' ' 
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one can write down the one-loop RG equations as 



oci = -ha?, (1) 
Yk = Y k (J2c ki ai -J2a M Yi), (2) 

i I 

where ■ = d/dt, t = log Mq UT /Q 2 and 
h = {33/5,1,-3}, 

c« = {13/15,3,16/3}, Cbi — {7/15, 3, 16/3}, c Ti = {9/5, 3, 0}, 
a u = {6,1,0}, a bi = {1,6,1}, a ri = {0,3,4}. 

The general solution to eqs.(|l],|D can be written as 



Y k °u k 
1 + a-kkY k ° Jo u k 



n = ; . \; n fc rt , (4) 



where the functions u k obey the integral system of equations 

E t E b E T 

Ut = 71 . ^rn rt. 777^5 U b = 7Z . „rn rt m/ c/ . . „m rt. 7777 ) U T 



(i + 6n°/o^) 1/6 ' (i + 6y t % t « t )V6(i + 4 y T o/ ^ T )V4' - (l + enUo^) 1 / 2 ' 

(5) 

and the functions E k are given by 

£?* = n(l + 6 i a?t) c «/ 6 '. (6) 

i=l 

Let us stress that eqs.(||,|) give the exact solution to eqs.(0,|]), while the u k s in eqs.(|5|), 
although solved formally in terms of the E k s and Y k 0, s as continued integrated fractions, 
should in practice be solved iteratively. Yet the important gain here is twofold: 

i) as shown in 0, the convergence of the successive iterations to the exact solution 
can be fully controlled analytically in terms of the initial values of the Yukawas, allowing 
in practice to obtain approximations at the level of the percent or less after one or two 
iterations and 

ii) the structure of the solutions is explicit enough to allow for exact statements about 
some regimes of the initial conditions, as we will see later on. Furthermore, these nice 
features will be naturally passed on to the solutions for the soft SUSY breaking parameters 
since the latter will be obtained from through the method of Ref.0. 



3 The Soft Terms and Grassmannian Expansion 

An important feature of the solution (|J[5p is that it is written in an analytic form with 
the initial conditions explicitly present. This allows one to get the same type of solution 
for all the soft terms in an iterative form. 
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Let us describe briefly the procedure. It has been shown[[8| that the soft terms which 
break supersymmetry can be introduced in a classical Lagrangian via the so-called spurion 
fields. This leads to the modification of the original couplings of a rigid theory, they 
become external spurion superfields depending on Grassmannian parameters 1 . In the 
MSSM it looks like 

a — > o.i = ai(l + mifj + rhifj + 2rni7fii7]f]) , (7) 
Y k ^Y k = Y k (l- A k T]- A fc r/ + (S fc + A k A k )rif}), (8) 

where are the gaugino masses, A k are the scalar triple couplings and are certain 
combinations of the soft masses 

= rh 2 Q3 + rh 2 U3 + m 2 H2 , E 6 = rh 2 Q3 + m 2 D3 + m 2 H1 , S r = m| 3 + m| 3 + m 2 Hl . 

Here rj = 6 2 and fj = 2 are the spurion fields depending on Grassmannian parameters 

e a ,e a ( a = 1,2). 

It has been proven in Ref. |J that the singular part of effective action, which determines 
the renormalization properties of any softly broken SUSY theory, is equal to that of an 
unbroken one in presence of external spurion superfields. This means that in order to 
calculate it in a softly broken case one just has to take the unbroken one, replace the 
couplings according to eqs.(|7]||) and expand over the Grassmannian parameters rj and fj. 

Moreover, as it has been demonstrated in JF], the same replacement can be done 
directly in RG equations in order to get the corresponding equations for the soft terms, 
or even in the solutions to these equations. In the last case one obtains the solutions to 
the RG equations for the soft terms. Below we demonstrate how this procedure works in 
case of the MSSM, when combined with the solutions (H,f|). 



4 Analytical Solution to RG Equations for the Soft 
Terms 

Let us now perform the substitution (|7j,^p in (0-|^) and expand over r] and fj. Then the 
linear term in r\ will give us the solution for m« and A k and the i]fj terms the ones for 
Efe. (For simplicity, we do not consider here CP-violating effects and take all the soft 
parameters to be real valued.) The resulting exact solutions look similar to those for the 
rigid couplings 

, c A° k /Y k ° + a kk Ju k e k 

A k = -e k A . , (1U) 

l/y fc u + a kk J u k 



1 The resulting formulae coincide with those of Ref. |I(J except for some minor difference in higher 
loops. Since we consider only one-loop RG equations we ignore this difference here. Similar results were 
obtained also in Ref. Jill 
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£ fc + A\ + 2e k A k - 



(Al) 2 /Y°-Xl/Y k ° + a kk fu k Z k 



(11) 



1/Y k ° + a kk Ju k 

where the new functions e k and £ k have been introduced which obey the iteration equations 



1 d£ 4 Ag J u b - J w 6 e fc 



£ t d V l/Y b ° + 6ju b ' 
1 dE 6 A® J Ut — J u t e t A® J u T J u T e 7 



E b d V l/Y t ° + 6Ju t l/Y T ° + 4ju 7 



1 dE T t ^A° b J u b - / u b e b 



E T drj l/Y b ° + 6ju b ' 
1 d 2 E t 2 1 dE t A° b fu b -fu b e b ? 



E t dr)dfj E t d V l/Y b ° + 6ju b 



' A° b j u b - ju b e b 

t i/y 6 ° + 6/ U6 , 



((Eg + {At) 2 ) ju b - 2A° b J u b e b + J u b ^ / + 6 J 



1 d 2 E b 1 dE b 
+ 2- 



A° t J u t - J u t e t A° T J u T - J u T e T 



u b 



l/Y t ° + 6ju t l/Y? + Aju 7 
( A° T Ju T - Ju T e T \ 2 



E b drjdrj E b drj 

( A° t J u t - J u t eX 
+ ^ l/Y t ° + 6Jut ) ' "\ l/Y? + 4fu T J 
( A\ J u t - J u t e t \ ( A° T J u T - J u T e T \ 
{ l/Y t ° + 6ju t J { l/Y? + 4fu T J 

- ((E? + (A° t ) 2 ) Ju t - 2A° t J utet + J «&) / (± + 6 

- ((£° + (A° T ) 2 ) Ju T - 2A° T J u T e T + J «^ T ) / (± + A J u T 

1 rf 2 E T 1 dE T A° b J u b - J u b e b ( A° b J u b - J u b e b \ 2 
E T dridrj E T drj l/Y b ° + QJu b { l/Y b ° + QJu b J 

-3 ((E° + (A°) 2 ) |« 6 - 2A° | u b e fe + j u^) / + 6 / ^ 



(12) 



where the variations of E k should be taken at r) = r\ = and are given by 



1 dE k 



E k drj 
1 d 2 E k 



E k drjdrj 



v,n=o 



r],r]=0 



i=l 



(3 \ 3 3 

c ki aim\ + 2t Y c fci ai(m°) 2 - t 2 ^ c ki bia 2 {m 0y2 
i=l / i=l i=l 



(13) 
(14) 



When solving eqs.(|5|) and (12) in the n-th iteration one has to substitute in the r.h.s. the 
(n — l)-th iterative solution for all the corresponding functions. 
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In the particular case where Yb = Y T = eqs. fl5r|P^) give an exact and well known 
solutions already in the first iteration. 

Let us finally note that upon inspection of the solutions (|9|-|T2|), the A^s and Sj's 
depend respectively linearly and quadratically on the initial conditions, as expected, and 
thus can be generally cast in the form: 

A, b ,At)= E aj(t)A?+ E h(t)ml (15) 

j=t,b,r fe=l,2,3 

*W(*) = E E d ij (t)A^A° j + E %(t)A°m° 

i,j=l,2,3 i,j=t,b,r i=t,b,T,j=l,2,3 

+ k t (t)Y? t + fc 6 (t)£g + fc r (t)E°, (16) 

where the various running coefficients dj, c^, dij, and ki are fully determined by our 
solutions and can be seen to depend exclusively on the initial conditions of the gauge and 
Yukawa couplings. In Sec. 6 we will evaluate these coefficients at the E.W. scale, using a 
truncation of the general solutions. 



5 Quasi-Fixed Points and the Independence of Initial 
Conditions 

The solutions (@-|5|, ^HH) en j°y the nice property of exhibiting the explicit dependence 
on initial conditions and one can trace this dependence all the way down to the final 
results. This is of special importance for the non-universal case since one can see which 
of the parameters is essential and which is washed out during the evolution. In particular 
the solution for the Yukawa couplings exhibit the fixed point behaviour when the initial 
values are large enough. More precisely, in the regime F t , Y b , Y r ° — ► oo with fixed finite 
ratios Y t °/Y b ° = r 1; Y b °/Y® = r 2 , it is legitimate to drop 1 in the denominators of eqs.(f|, 0) 
(see appendix A for a proof) in which case the exact Yukawa solutions go to the so-called 
IR quasi-fixed points (IRQFP) defined by 



u 



FP 



yr = -^P (17) 



FP 

k 



with 



fp _ Et FP „.FP 



u b ( r..FP\M(\f C..FPM/4' u t /r„.FP\l/9 v xc 



t ~ {Ju FPy/^ 6 {Iu FPy /e{Iu FPy/*i r - {Ju FPy/2 

extending the IRQFP |l^] to three Yukawa couplings. What is worth stressing here is 
that both the dependence on the initial condition for each Yukawa as well as the effect of 
Yukawa non- unification, 7*i,r2 have completely dropped out of the runnings. ( Note that 
in practice this regime is already obtained if Y fc ° > aQ UT , assuming here for simplicity 
the unification of the gauge couplings a® = a° = a° = o>q UT .) The fact that the ratios 
ri,r 2 drop out implies the validity of the described properties in any tan/? regime. 
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This in turn leads to the IRQFPs for the soft terms. Disappearance of Y k ° in the FP 
solution naturally leads to the disappearance of A k and E° in the soft term fixed points. 
To see this one can either go to the limit of large Y k in eqs. (|9|-|T2]) or directly perform the 
Grassmannian expansion of the FP solutions ([17|JT8| ). One gets 

s<^r 

s< 



A F k P = -e! P + J -^-Fh, (19) 



V FP C FP , f A FP\2 , FP A FP J a k Sfc 

^k — Sfc + {^k ) ' ze k A k r I FP 



Jut 



r FPcFP ffii FP p FP \ 
rFP _ J u k Z,k _ ( p FP\2 , f J u k e k / 9f) x 



with 



e FP 



FP 



e FP 



1 dE t 






E t d V 


e j< p ' 




1 dE b 


i surer 




E b dr] 


e j u r 




1 dE T 


i surer 




E T drj 


2 jur ' 





[21) 



i ^ , 7 f jur ery 1 1 ^/< p e f p i/«rer 



^ E t dridfj 36\ Ju FP J 3E t dri J u pp 6 /wf p 

1 7 //n^^ 



Eb drjdfj 36 

1 d£> 



1 o? -£/ T 
^ r ~ E~dr]df] 4 




One can see from eqs. ( |T9| , pO|) that the constant terms in and do not contribute to 
Afc and and can be dropped from eqs. fl2l"|) . Thus, all the dependence on the initial 
conditions Y k , A k and TP k disappears from the fixed point solutions. The only dependence 
left is on the gaugino masses. This is a general screening property valid for the exact 
solution as well as for any of its approximate (truncated) forms. 

In view of the above screening properties of the initial conditions at the quasi-fixed 
point, one should recall the existing connection between the true IR attractive fixed point 



of the Yukawa couplings and that of the soft parameters ||13||. What we have established 



can be seen as an extension of such connections to the transient regime of quasi-fixed 
point at the one-loop level. It is also worth stressing that the above properties are valid 
for any renormalization scale and are thus operative despite the uncertainty in the choice 
of the physical scales. 
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6 Numerical Analysis. The Role of Non-universality 



We perform now a numerical study of our solutions in the first and next iterations in order 
to demonstrate the convergence of the procedure and to show the role of certain initial 
conditions. In particular we will show that assuming that the soft terms are of the same 
order of magnitude at the GUT scale the only essential parameter is the gluino mass. All 
the other soft terms are suppressed either due to the fixed point behaviour mentioned 
above, or due to the smallness of ot\ and a 2 compared to a% at low energy. 



6.1 The first iteration 

To have a solution in the first iteration one takes eqs.([|,[l^) where in the r.h.s. the 
functions u k , e k and are taken in the 0-th iteration. They are 

1 dE k 1 d 2 E k . . 

u k = E k , e k = — — , £ fc = — — , 22 
hi k arj ti k dr/drj 



where the variations of E k are given by eqs.(I3| 
For comparison we consider two particular cases: 

i)F, o = 5a , n° = n° = o, 

II) Y t ° = Y b ° = F T ° = (14- 10)« , ( with a « 0.00329). 

The first case corresponds to the so-called low tan (3 regime, where the first iteration 
is already exact, and the second case refers to the high tan/3 regime for an 5O(10)-like 
initial conditions. 

In the first case the only essential soft terms besides the gaugino masses are A t and 
E t . All the soft masses of the third generation are expressed through £ f in a simple 
way[|IJ]. A t and H t have a dimension of a mass and mass squared, respectively, and may 



be expanded over the initial conditions in the following fashion (eqs.([T5|, |T6|) , 

A t = 0.0359 A t0 - 0.0167 m w - 0.1598 m 20 - 1.6284 m 30 , (23) 
£ t = -0.0346 A 2 t0 + 0.0044 A to m w + 0.0270 A t0 m 20 + 0.1225 A tQ m 30 

+0.0184 m 2 10 - 0.0057 m 10 m 2 o + 0.2566 m 2 20 

—0.0335 mi m 30 — 0.2735 77120^30 + 6.3695 m\ 

+0.0359 (m° H2 ) 2 + 0.0359 (mj 3 ) 2 + 0.0359 (m^) 2 , (24) 

where the numbers are calculated for t = log Mq UT /M^ w 66. 

One can see from eqs.(^,^) that the prevailing term is that of m\. A® and E° 
decouple due to the fixed point behaviour as explained above and the contribution of m\ 
and 777.2 ls sma U due to smallness of a.\ and a 2 compared to 0:3 at t — 66. 

The second case looks similar. We first consider the triple couplings A k . In this case 
one has a set of initial values {A^m^}. In fig.l we show the variation of the coefficients 
di, b k of eq.flTo]) for A t , A^ and A T as functions of Y k in the interval Y k = a + 10«o 

One can clearly see that the coefficients of A k are small and have a fast decrease with 
increasing Y k °. The coefficients of m" quickly saturate and approach their asymptotic 
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values with the hierarchy 1 : 10 : 100 for m\,m 2 and m". The effect is less pronounced 
for A T due to the absence of the SU(3) coupling in the lepton sector. 

Now we come to £*.. We have chosen the intermediate value of Yj? = 5ao where the 
effective fixed point is practically already reached and calculate the coefficients at t — 66 
as in eqs. (|23|j24] ). One has 

£ t = -0.0497 A 2 t0 + 0.0076 A 2 m + 0.0142 A t0 A b0 + 0.0057 A to m w 

-0.0018 A b0 m l0 + 0.0333 A t0 m 20 - 0.0114 A b0 m 20 + 0.1509 A t0 m 30 
-0.0516 A b0 m 30 + 0.0198 m\ Q + 0.2509 m 2 + 6.3299 m 2 30 
-0.0057 m 10 m 20 - 0.0336 m 10 m 30 - 0.26 69 m 20 m 30 - 0.0252 {m° Da ) 2 
-0.0252 (m° Hi ) 2 + 0.0525 (m° H2 ) 2 + 0.0273 (m° Q3 ) 2 + 0.0525 (m^ 3 ) 2 , (25) 



S fe = +0.0079 A 2 t0 - 0.0717 A 2 b0 + 0.0058 A 2 r0 + 0.0200 A m A m 

-0.0062 A t0 A r0 + 0.0262 A b0 A r0 - 0.0005 A to m w + 0.0013 A bo m w 
+0.0004 A TO m 10 - 0.0119 A t0 m 20 + 0.0395 A bQ m 2Q - 0.0121 A T0 m 20 
-0.0618 A to m 3 o + 0.2008 A bo m 30 - 0.0965 A T0 m 30 - 0.0052 mf 
+0.2359 m 2 20 + 6.8165 m^ - 0.0027 m w m 20 - 0.0037 m w m 30 
-0.2498 m 20 m 30 + 0.0778 (m° Ds ) 2 - 0.0502 (m^ 3 ) 2 + 0.0276 (m° Hi ) 2 
-0.0324 (m° H2 ) 2 - 0.0502 (m?J 2 + 0.0454 (mjj 2 - 0.0325 (m Ua ) 2 , (26) 

S T = —0.0010 Al — 0.1947 A\ + 0.1397 A M ^4 T o — 0.0195 A bo m w 

+0.0339 A TO m 10 - 0.0409 A b0 m 20 + 0.0718 v4 r0 m 20 + 0.0894 A b0 m 30 
-0.0842 A rQ m 3Q + 0.0993m 2 + 0.3527 m 2 - 2.8162 m 2 
-0.0133 mi m 2 o + 0.0116 m 10 m 3 o - 0.0687 m 20 m m - 0.2177 (m^) 2 
+0.2649 (m Ea ) 2 + 0.0473 (m^) 2 + 0.2649 (m$J 2 - 0.2177 (my 2 . (27) 

One can again see how the coefficients of the initial values of A® and E° almost vanish 
and the prevailing one is that of (m|}) 2 . The next-to-leading ones are those of (m 2 ) 2 and 
m 2 m3 being however almost 30 times smaller. This is true for both E° and S° but is less 
manifest for E°. We note here that a soft gaugino mass hierarchy like the one predicted 
by anomaly- mediated susy breaking, : m 2 : ml = 3 : 0.3 : 1 @, enforces even more the 
insensitivity of the running Aj's and Ej's to the non- universality of the gaugino sector. 



6.2 The next iterations 

To demonstrate the validity of the iterative procedure and reliability of the the first 
iteration we consider the effect of the next ones on the above mentioned coefficients. We 
have performed the numerical integration up to the 6-th iteration and have observed fast 
convergence of the coefficients to their exact values. To show the numbers we have chosen 
the leading coefficients of tuq 3 in A^ and rn^ 3 in £*.. In case = Y b ° = Y® = 5ao and 
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«o = 0.00329 the results are the following: 



Iteration 


A t 


A b 




A T 




Eft S T 




-1.6127 


-1.7584 





6871 


6.3299 


6.8166 -2.8162 


<2nd 


-1.6161 


-1.7330 





5037 


6.3270 


6.6822 -2.0989 


3 rd 


-1.6125 


-1.7372 





5526 


6.3192 


6.7069 -2.2937 




-1.6133 


-1.7375 





5440 


6.3213 


6.7054 -2.2588 


5 th 


-1.6131 


-1.7373 





5456 


6.3206 


6.7053 -2.2660 


Qth 


-1.6131 


-1.7374 





5454 


6.3207 


6.7054 -2.2649 



One can see explicitly the fast convergence of the iterations. As expected it is worse 
for A T and E r , so in this case one has to take few more iterations. We present the 
general arguments for the convergence of iterations for the soft terms in appendix B. The 
advantage of this solution is that one can improve the precision taking further iterations 
and in principle can achieve any desirable accuracy. Typically one has an accuracy of 
a few percent after 2-3 iterations. This is in contrast with the approximate solutions 



presented in Ref.jnj] which give simple explicit expressions but without improvement. 
Taking the sixth iteration in eqs. (|5|JT^) expressions for the soft terms now look like 

A t = 0.0558 A t0 - 0.0294 A m + 0.0080 A r0 - 0.0186 m 10 - 0.1586 m 20 - 1.6131 m 30 , 
A b = -0.0341 A t0 + 0.0984 A b0 - 0.0450 A r0 - 0.0014 m 10 - 0.1583 m 20 - 1.7374 m 30 , 
A T = 0.0394 Ato - 0.2221 A b0 + 0.2871 A rQ - 0.0825 m 10 - 0.2344 m 20 + 0.5454 m 30 , 



£ t = -0.0487 A 2 + 0.0068 A 2 m + 0.0014 A 2 t0 + 0.0130 A m A b0 

-0.0017 A t0 A r0 - 0.0017 A b0 A r0 + 0.0058 A to m w - 0.0015 A bo m w 
-0.0001 A r0 m 10 + 0.0335 A to m 2 o - 0.0097 A w m 20 + 0.0005 A T0 m 20 
+0.1514 A t0 m 30 - 0.0460 A w m 30 + 0.0070 A rQ m 30 + 0.0211 m 2 w 
+0.2547 m 2 20 + 6.3207 m 2 30 - 0.0057 m 10 m 2 o - 0.0340 m 10 m 30 
-0.2720 m 20 m 30 - 0.0294 (m^) 2 - 0.0214 (m° Hi ) 2 + 0.0558 (m^) 2 
+0.0264 (m° Qs ) 2 + 0.0558 (m^ 3 ) 2 + 0.0080 (m^) 2 + 0.0080 (m^ 3 ) 2 , 

E 6 = +0.0067 A 2 - 0.0736 A 2 b0 - 0.0035 A 2 r0 + 0.0200 A t0 A w 

-0.0064 A t0 A T0 + 0.0302 A b0 A r0 - 0.0003 A to m 10 - 0.0000 A bo m 10 
+0.0021 A r0 m w - 0.0110 A t0 m 20 + 0.0396 A b0 m 20 - 0.0076 A r0 m 20 
-0.0613 A t0 m 30 + 0.2356 A b0 m 30 - 0.0945 A r0 m 30 - 0.0013 m 2 10 
+0.2532 m 2 + 6.7053 m 2 - 0.0030 m 10 m 20 - 0.0053 m 10 m 30 
-0.2578 m 20 m 30 + 0.0984 (m^) 2 + 0.0534 (m° Hi ) 2 - 0.0341 (m^) 2 
+0.0642 (m° ft ) 2 - 0.0341 (m^) 2 - 0.0450 (m£ 3 ) 2 - 0.0450 (m^ 3 ) 2 , 

S r = +0.0009 A 2 t0 - 0.0106 A 2 m - 0.1862 A 2 + 0.00 77 A t0 A b0 



10 



-0.0146 A t0 A T0 + 0.1334 A b0 A T0 + 0.0021 A to m 10 - 0.0190 A w m 10 
+0.0350 A r0 m 10 + 0.0007 A tQ m 2 o - 0.0350 A b0 m 20 + 0.0743 A T0 m 20 
-0.0397 A t0 m 30 + 0.1288 A b0 m 3Q - 0.1090 A r0 m 30 + 0.1029 m 2 
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+0.3907 m 2 20 - 2.2649 m 2 - 0.0140 m 10 m 20 + 0.0160 m 10 m 3 o 
-0.0504 m 2Q m 30 - 0.2221 (m^) 2 + 0.0650 (m° Hi ) 2 + 0.0394 (m° H2 ) 2 
-0.1827 (m° ft ) 2 + 0.0394 (m^) 2 + 0.2871 (m£ 3 ) 2 + 0.2871 (m° E3 ) 2 



to be compared with eqs.(25,26,2^). These numbers can be now used for the calculation of 
the soft parameters and masses using eqs.(^) with arbitrary initial conditions for the soft 
terms. Being calculated for Y t ° = Y b ° = Y® = 5a and a — 0.00329 they are practically 
independent of the initial values of Yukawa couplings provided the latter are big enough. 



7 Towards the Physical Masses 

The values of E& completely define those of the soft masses for squarks, sleptons and 
Higgses due to linear relations which follow from the RG equations |nj] and read, after 
relaxing the universality assumption^, 



o 2 128/3 + 87/ 2 - , 17(S t -Eg)+20(E b -Eg)-5(E T -E?) 
m Q3 = (m Q3 ) + — + — , (28) 

~ 2 ,-o x 2j 144/3 - 108/2 + 144/5/! 
m U3 = («% 3 ) + — + 



- 2 ,~o^ 2l 144/ 3 - 108/ 2 + 144/5A , 42(E t - E°) - 8(E b - Eg) + 2(E T - E^ 



122 122 

,~o , 2 , 112/3-84/2 + 112/5/! -8(E f - Eg + 48(E b - Eg) - 12(E r - Ej) 
m m = m m H 

D3 V D3J -r ^ 122 

240/ 3 - 3/ 2 - 57/5/i , -9(S, - E t °) + 54(E 6 - Eg) + 17(E T - E°) 



m v H1J 122 122 

a -272/3 + 21/2 - 89/5A 63(S f - Eg) - 12(E b - Eg) + 3(E r - Eg 

m H2 - [m H2 ) + — + — 

-2 ,-0 , 2j 8O/3 + 123/ 2 - 103/5A 3(E f -Eg)-18(E b -Sg) + 35(E r -E°) 

m r , = m r o H , 

L3 V L3; -r t 122 

_ 2 o , 2 I6O/3 - 120/2 + 32/i , 6(E t -Eg)-36(E b -Sg) + 70(E r -E 

™£3 = ( m i?3) + ~ + 



0> 



where 



122 122 



hi \ (l + bia t) 2 J 

At this level one can already make rough qualitative statements about the physical 
scalar masses. We note first that, as can be seen from the above equations, the sensitivity 
to the initial conditions reappears partly in the running of the soft scalar masses, even in 
the vicinity of the IRQFP. However, this dependence remains confined in the initial values 



2 Note that even though a trace term "Tr(Y/ l j, perc 7 iar . ge m 2 )" is generically present in the RGE in the 
non- universal case, it cancels out in eq.(p8|) 
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of the soft masses themselves in a universal (scale independent) form, and in the initial 
conditions of the gaugino soft masses through the /j's and the E's. [The dependence on 
the initial values of Yukawa couplings as well as on the v4°'s and the E°'s, that could come 
from the running, remain completely screened.] The ratios giving the universal sensitivity 
of the running soft scalar masses to the soft scalar masses initial conditions is as follows: 







« 3 ) 2 : 


(rh° D3 ) 2 : 


Km) 2 : 


Km) 2 : 


K 3 ) 2 : 


K E3 ) 2 


(™Qs) 2 
(m U3 ) 2 
(m D3 ) 2 
(mHi) 2 
{m H2 ) 2 
Ks) 2 
{m E3 ) 2 


17 


-3.4 


-4 


-3 


-3.4 


1 


1 


-17 


40 


4 


3 


-21 


-1 


-1 


-5 


1 


9.25 


-4.5 


1 


1.5 


1.5 


-5 


1 


-6 


5.67 


1 


-1.89 


-1.89 


-17 


-21 


4 


3 


19.67 


-1 


-1 


5 


-1 


6 


-5.67 


-1 


29 


-11.67 


5 


-1 


6 


-5.67 


-1 


-11.67 


8.67 



These numbers are renormalization scale independent and give the trend of the relative 
sensitivity in the vicinity of the IRQFP. 

On the other hand, the dependence on the initial soft gaugino masses is renormaliza- 
tion scale dependent. At the electroweak scale (t ~ 66), one finds that in the soft masses 
of the third squark generation and of the Higgs doublets the sensitivity to (m®) 2 remains 
leading (by a factor of 15 to 25) as compared to (m®) 2 . In contrast, a large cancellation 
occurs for the sleptons, leading to comparable sensitivities to (m\) 2 and (m^) 2 in rfiE3, 
and even a bigger sensitivity to (m") 2 (by a factor of 4) in m E3 . 



To go further to the physical scalar masses, one has to consider the behaviour of the 
fi parameter which enters the mixing of the left and right states. The running of this 
parameter has the simple form //(t) ~ fioexplf^a — Y)] where a, Y are generic gauge and 
Yukawa couplings. Thus, here too, the initial conditions for the Yukawas are screened 
near the IRQFP in the evolution of //, the Aj's and Ej's being absent anyway. However, 
when the electroweak symmetry breaking (EWSB) is required to take place radiatively, 
the /1 parameter becomes, as usual, correlated to the other parameters of the MSSM at 
the electroweak scale. To be specific, in the leading one-loop top/stop-bottom/sbottom 
approximation to the EWSB conditions, the sensitivity of fi to the initial conditions will 
come basically from the soft scalar masses of the Higgs doublets and scalar partners of 



the third quark generation |T5|. As stated before, the latter dependence is dominated, 
on one hand by the initial conditions of the soft scalar masses, in a well determined scale 
independent way, and on the other hand by the (scale-dependent) m 3 contributions. The 
same dependence pattern is then taken over to the physical scalar masses. A further in- 
clusion of the scalar r contributions to the EWSB conditions will basically not affect this 
dependence pattern. Indeed, although fh 2 L3 and fh 2 E3 have comparable sensitivity between 
and mg at the electroweak scale, they are less sensitive to this sector altogether than 
to the squark soft masses. All in all our analytical results allow to draw at this stage a 
qualitative sensitivity hierarchy for the physical scalar masses: 
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- basically no sensitivity to Yukawa couplings initial conditions (whether unified or 
not), or A® initial conditions (whether universal or not), 

- important sensitivity to initial conditions of the soft gaugino masses, however basi- 
cally only through 777.3, J - a we& k sensitivity to non-universality of this sector, 

- important sensitivity to initial conditions of the soft scalar masses, however through 
a universal scale independent pattern. 

8 Conclusion 

In the present paper we have obtained general analytical forms for the solutions of the 
one-loop renormalization group equations in the top/bottom/r sector of the MSSM. These 
solutions are valid for any value of tan (3 as well as any non-universal initial conditions for 
the soft SUSY breaking parameters and non-unification of the Yukawa couplings. They 
allow a general study of the evolution of the various parameters of the MSSM and to trace 
back, sector- wise, the sensitivity to initial conditions of the Yukawa couplings and the soft 
susy breaking parameters. We have established analytically a generic screening of non- 
universality, in the vicinity of the infrared quasi fixed points. In practice, this property 
gives the general trend of the behaviour, despite the large number of free parameters, 
and even when one is not very close to such a quasi fixed point. This shows that non- 
universality of the A parameters and gaugino soft masses, as well as Yukawa unification 
conditions, would basically have no influence on the squark and Higgs spectra. The 
main input from the gaugino sector comes from the soft gluino mass contribution (which 
dominates by far the other two), i.e. insensitive to non-universality conditions of this 
sector. The only substantial sensitivity to non-universality is associated to the initial 
conditions of the scalar soft masses, but is renormalization scale independent and well 
defined. A similar pattern holds for the sleptons, apart from the fact that now the 
contribution of the wino soft mass becomes comparable to that of the gluino, yet the 
overall sensitivity to the gaugino sector is much smaller than in the case of the squarks. 

Detailed illustrations of the physical spectrum, including the lightest Higgs, will be 
given in a subsequent study. 
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Appendix A: Screening of initial conditions and of 
non-universality at the IRQFP 

We give here a short proof that in the regime Y®, Y b °, Y® — > oo with fixed finite ratios 
Y®/Y b ° = ri,Y®/Y® = r 2 the Yukawas become insensitive both to Y k ° and r 1; r 2 . This 
result would be immediate from eq.(|]) by dropping 1 in the denominator were it not for 
the fact that the u k s have also a non-trivial dependence on the large initial conditions. 
Let us rewrite eqs.([D in the form 



~(n) 



~{n) 



;i + 6n°/4 n " 1) ) 1/6 

Eh 



(1 + 6Y b ° J ui^y/^l + AY b ° J ^ n " 1) ) 1/4 



4^ = > (A-l; 

(i + en° / uty /2 



with 



ut ] =E k , (k = t,b,r), 

where the twiddled quantities are obtained from the non twiddled ones by proper rescaling 
with r\ or r 2 , and we indicate explicitly the order of iteration. It is now easy to show 
inductively that if at the n th iteration 



~,( n )\FP 

for Y b ° oo, with < p k < 1 , (A.2) 



~(n) {Uk V ' f_ vO 
k (Y b °)P* 



where (u k nS> ) FP is Y b independent but r 1; r 2 dependent, then the same is true at the 
[n + l) th iteration. Furthermore, since ( |A.2|) is obviously true for n = 1 as can be easily 



seen from (|A.1|) we conclude that the exact u k s behave also like 



(Uk ) 



FP 



u k ~ TTToWk with < Pk < !■ 



This means that the l's can be legitimately dropped both in eq.@ and eq.(pj). The 
complete cancellation of Y b , r\ and r 2 in the final result is then obvious, leading to eqs. (|T7| , 

m). 



Appendix B: Convergence of iterations for the soft 
terms 

In this appendix we prove that the convergence of the e^s and is automatic once that 
of the u'iS is achieved. In particular this means that a controllable behaviour is expected 
whatever the initial conditions for the soft parameters may be. 
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Let us define 



e t {t) 
S{t) = | e b (t) 
e T (t) 

-U b (h;t) 

U(h;t) = | -U t {t x ;t) -U T {h;t) 
-3U b (ti;t) 



where 



ZYjfti; i) = — ; — "' v " ±y — i = t,b,r 
1 ' 1/1? + 0^' 



and at = a b = 6, a T = 4, 



(B.l) 



(B.2) 



(B.3) 



C(t) 













(B.4) 

with C(O) = 0. The system of integral equations for the e^'s can then be written in the 
matrix form 

S(t) = C(t) + [ U{tx;t)E(tx)dbi (B.5) 
Jo 

To prove the convergence of S(t) we the define the mapping £ — > £': 

£'(t) = C(t) + f U(h; fiettjdh (B.6) 
Jo 

and the norm | . | through 



|| M(t) ||= sup {max | |} 



0<t<T 



(B.7) 



for any matrix M in a given evolution interval [0, T]. One then has the inequality 

I f{U£)ij\<(£f\ U ^\)\\ £ \\ (B.8) 
Jo k Jo 

valid for any On the other hand, one has from eqs.( B.2| , |B.3| ) 



ik — S 



I u b < 1 

i/y 6 °+6 J' Ufc - 6 



(* = 1) 

.,._,.< S ( < = 2 ) 

3 l/Y°+6fu b -5 (* = 3 ) 



l/y t °+6 Jut l/K T °+4 / u T 



(B.9) 
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Combining the above inequalities ( |B.8| , |B.9| ) with eq. ( |B.6| ) one obtains 



£■[ — £'2 II < 2 II ^ l ~ ^ 2 



(B.10) 



that is, the mapping ( |B.6[ ) is a contraction, the solution to eq.( |B.5| ) is unique and ap- 
proximated at worse with an error of l/2 n , after n iterations. Actually, the situation is 
much better than given by this upper bound error, as one can see from the numerical 
illustrations of section 6.2. Finally, we note that the rational is exactly the same for the 
convergence of the £'s. Indeed, apart from a different definition for C(t), the £'s satisfy a 
matrix equation similar to ( p.5| ) with the same U as the one given in (P-2| ). 
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Figure 1: Initial value contributions of the various soft SUSY breaking parameters to 
the running A t , A b et A T at the EW scale, as a function of a common initial value for the 
three Yukawa couplings 
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